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Abstract. Supermassive black holes are a key element in our understanding of how galaxies 
form. Most of the progress in this very active field of research is based on just ~ 30 determi- 
nations of black hole mass, accumulated over the past decade. We illustrate how integral-field 
spectroscopy, and in particular our OASIS modeling effort, can help improve the current situation. 
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1. Supermassive black holes and galaxy evolution 

Fifteen years ago the existence of supermassive black holes (BHs) in galaxy nuclei was 
considered an interesting possibility which had to be demonstrated. Nowadays BHs are 
regarded as the crucial ingredient for our understanding of how galaxies form. Key to this 
paradigm shift was the launch in 1990 of the Hubble Space Telescope (HST). It all started 
with the realisation that the mass of the BH is correlated to other global characteristics 
of the host galaxy as a whole. Initially a correlation Mbh — L was found between the mass 
of the BH and the luminosity of the host-galaxy stellar spheroid ( Kormendv fc Richstond 
1995t iMagorrian et ai1ll998th In 1997 the installation of the STIS long-slit spectrograph 
on HST allowed the spatially-resolved kinematical observations to probe inside the radius 
of the subarcsecond BH sphere of influence Rbh = GM^n/ ' 2 m nearby galaxies (a 
being the velocity dispersion of the galaxy stars). The increased accuracy in the Mbh 
determinations contributed to the discovery of the much tighter Mbh — o correlation 
(|Gebhardt et alj|200fi iFerrarese fc Merrittll2000l ). 

Simil ar correlations were f ound between the Mb h and respectively the galaxy concen - 
tration (| Graham et al. 200lh. the dark-halo mass (lFerrare"sel 12002 ; Pizzella et al. 2005 ). 
the galaxy mass ( Marcon i & Hunt 2003t lHaring k, Rixll2004 ) and the stars gravitational 
binding energy (|Aller fc R,ichstone!l2007t h The existence of these correlations is broadly 
consistent with a scenario in which the BH regulates the galaxy formation, during the 
hierarchical galaxy merging, by shu tting off the conversion of gas into stars vi a a feedback 
mechanism due to its powerful jet ( Silk &: Reesl 19981 Di Matteo et al. 2005 ). 
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Figure 1. Updated Mbh — cr relation. Open squares are values from the literature; open circles 
are literature val ues, plotted against th e luminosity- weighted a within lR e determined from 
SAUR0N data (see ICappellari et al J 1200 6); filled circles are our new Mbh OASIS de t ermin ations, 
against the SAURDN a. For reference, the solid li ne is the relation of lTremaine et aD ((20021 ) . while 
the dotted one is from iFerrarese fe Fordl (|2005l ) . 



2. Observational evidences and current limitations 



Our understanding of the role of BHs in galaxy evolution is however far from complete. 
In fact the current observables do not uniquely constrain the models, which depends on 
a number of assumptions. This is due in part to the relatively small number of secure BH 
measurements: in a deca de of high-resolution ob servations and models only ~ 30 values 
have been obtained (see IFerrarese fc Fordl 120051 for a recent review) . It is remarkable 
that so much progress was based on an extrapolation of so few BH measurements, and 
for a biased galaxy sample, to the entire galaxy population! 

An additional complication comes from the fact that the above correlations provide a 
rather indirect test for the models. A complementary and more direct approach to test 
the BH formation paradigm consists of looking in nearby galaxy centres for the signatures 
of the joint formation of the BH and the galaxy spheroid. Simulations of galaxy mergers 
show in fact that, when two galaxies with a BH merge, the distribution of the stellar 
orbits in the resulting remnant, after the two BHs coalesce, is significantly different from 
the one of the progenitor galaxies. This is due to the ejection of stars passing, along 
radial orbits, close to the resulting BH binary. Two observables signatures are expected: 
(i) the density profile should flatten inside the core radius Rc , which is much larger than 
-Rbh and (ii) the orbital distribution should be biased toward s tangential orbits inside 
R c (|Quinlan fc Hernauist|[l997HMilosavlievic fc Merrittll200ll ). 

Evidence for the formation of the density core , and its expected relation to the BH 



mass , was found from photometric observations (|Faber et alj|1997t iMilosavlievic et al 



2002). The detection of the orbital signature is more co mplicated as it requires integra l- 
field spectroscopic observations at high resolution fe.g. ICappellari fc McDermidll2005l ). 
The orbital distribution, in a non-spherical stationary system, is in fact a function of the 
three isolating integrals of motion and requires at least a three-dimensional observable 
quantity to be constrained. Nonetheless a first attempt at deriving the nuclear stellar 
anisotropy, for a samp le of 12 galaxies, was done using long-slit STIS spectroscopy by 
Gebhardt et al.l ( 20031) . They found a tendency for the orbits of the most massive objects 
to be tangentially biased. However this appears to be true only well inside -Rbh and not 
up to Rc as predicted by t he simulations. Si milar result s were recently f o und f or another 
carefully studied galaxy by Houghton et al. ( 20061 ) and Gebhardt et al. ( 2007 ). 
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Figure 2. Anisotropy profiles from the dynamical models. In standard spherical coordi- 
nates or is the second moment of the velocity distribution along the radial direction, while 
erf. = (a'g +<rl) /2. Different lines represent measurements at equally-spaced radial sectors in the 
galaxy meridional plane, from the equatorial plane to the symmetry axis. The vertical dashed 
and solid thick line indicate the position of Rbh and Rc respectively. The values of Rc an d the 
logarithmic nuclear slope 7' in the surface brightness were taken from lLauer et all l|2007l ). 



3. The role of OASIS integral-field spectroscopy 

Ground-based and hig-S/N integral-field observations of the stellar kinematics can 
overcome the limitations discussed in the previous section, and can be obtained for large 
galaxy samples using large telescope mirrors. By tightly constraining the orbital struc- 
ture, dynamical models fitted to integral-field observations (i) can measure BH masses 
with accuracy comparable or better than that ob tained with HST spectroscopy, even 
when i?BH is not well resolved ( Shapiro et al. 20061 ): (ii) do not suffer from the degener- 
acy in the recovery of the orbital distribution. 

For this we observed with t he OASIS integral-field spectrograph a sample of 28 el- 
liptical and lenticul ar galaxies dMcDermid et ah 2006[ ). The galaxies were selected from 
the SAURON sample (|de Zeeuw et al.ll2002l ). for which the needed l arge-scale integral-fiel d 



kinematics is also available up to about one half-light radius R e (|Emsellem et a.l.M2004f ) . 
The OASIS observations complement the SAURON ones by providing an order of magni- 
tude increase in the pixels density and a factor of two improvement in the median seeing, 
resulting in subarcsecond resolution. 

Here we report some results of the stellar dynamical models for an initial set of eight 
galaxies from the OA SIS sample. We constructed the models using our axis ymmetric im- 
plem entation ( Cappellari et al.ll2006f) of the orbital-superposition method ( Schwarzschildl 
1979f) . and we combined the SAURON and OASIS kinematics as in lShapiro et all ( 2006t ). 
We find that the Mbh is recovered with a median formal error of 30% (at the 3rj confi- 
dence level, for one degree of freedom, i.e. A\ 2 = 9). This accuracy is similarly to that 
of the previous HST determinations using long-slit spectroscopy. Moreover our integral- 
field kinematics allows for a more rigorous and robust determination of the luminosity- 
weighted second velocity moment a inside R e , which is used in the Mbh — o diagram. Our 
new measurements appear to follow within the scatter of the previous values (Fig.[T]). 

The integral-field spectroscopy allows for a robust recovery of the nuclear orbital distri- 
bution. Our first results do not show any evidence for the predicted connection between 
the existence and strength of the density core, and the presence of tangentially biased 
orbits within Rc (Fig. [3J) . Galaxies with a flat core do not necessarily show an increase 
of the tangential orbits inside Rc- Tangential anisotropy is generally found only near 
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the smaller -Rbh, at the limit of our spatial resolution. Our recovered anisotropy is more 
robust and appears less noisy than that derived from previous models based on long-slit 
observations, but our finding is not inconsistent with them. 

A statistically more significant sample, with models based on integral-field observa- 
tions, is needed to confirm the apparent discrepancy between the observations and the 
predictions of the binary-BH core-scouring mechanism. Additional accurate BH determi- 
nations are needed to explore the dependence, of the various correlations involving BHs, 
on other galaxy global parameters. Our OASIS modeling effort provides a step in this 
direction. However higher-resolutio n integral- field obser vations, as can be obtained with 
the adaptive optics technique (e.g. iNowak et al.l 12007). are needed to study less mas- 
sive galaxies, and to probe deeper inside Rbh ■ At the same time new and more realistic 
N-body simulations, in a cosmological context which includes the hierarchical merging 
process, are required to be compared with the observations. This will ultimately confirm 
our current understanding of the role of BHs in galaxy formation or require a revision of 
the picture we constructed in the past decade. 



References 

Aller, M. C, & Richstone, D. O. 2007, ApJ, 665, 120 

Cappellari, M., & McDermid, R. M., 2005, Class. Quantum Grav., 22, 347 
Cappellari, M., et al., 2006, MNRAS, 366, 1126 
de Zeeuw, P. T., et al. 2002, MNRAS, 329, 513 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604 
Emsellem, E., et al. 2004, MNRAS, 352, 721 
Faber S. M., et al., 1997, AJ, 114, 1771 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Ferrarese, L. 2002, ApJ, 578, 90 

Ferrarese, L., & Ford, H. 2005, Space Science Reviews, 116, 523 
Gebhardt, K., et al. 2000, ApJ, 539, L13 
Gebhardt, K., et al. 2003, ApJ, 583, 92 

Gebhardt, K., et al. 2007, ApJ, in press l|arXiv:0709.0585t ri [astro-ph]) 
Graham, A. W., Erwin, P., Caon, N., & Trujillo, I. 2001, ApJ, 563, Lll 
Haring, N., & Rix, H.-W. 2004, ApJ, 604, L89 

Houghton, R. C. W., Magorrian, J., Sarzi, M., Thatte, N., Davies, R. L., & Krajnovic, D. 2006, 

MNRAS, 367, 2 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Lauer, T. R., et al. 2007, ApJ, 664, 226 
Magorrian, J., et al. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 
McDermid, R. M., et al. 2006, MNRAS, 373, 906 
Milosavljevic, M., & Merritt, D. 2001, ApJ, 563, 34 

Milosavljevic, M., Merritt, D., Rest, A., & van den Bosch, F. C. 2002, MNRAS, 331, L51 
Nowak, N., Saglia, R. P., Thomas, J., Bender, R., Pannella, M., Gebhardt, K., & Davies, R. I. 

2007, MNRAS, 379, 909 
Pizzella, A., Corsini, E. M., Dalla Bonta, E., Sarzi, M., Coccato, L., & Bertola, F. 2005, ApJ, 

631, 785 

Quinlan, G. D., & Hernquist, L. 1997, New Astronomy, 2, 533 
Schwarzschild, M. 1979, ApJ, 232, 236 

Shapiro, K. L., Cappellari, M., de Zeeuw, P. T., McDermid, R. M., Gebhardt, K., van den Bosch, 

R. C. E., Statler, T. S., 2006, MNRAS, 370, 559 
Silk, J., & Rees, M. J. 1998, A&A, 331, LI 
Tremaine, S., et al. 2002, ApJ, 574, 740 



